THE ALGEBRA OF FACTORIAL POLYNOMIALS 

D. BABUSCI, G. DATTOLI, AND M. CARPANESE 

Abstract. We discuss the formal aspects of the factorial polynomials and of the associated series. 
We develop the theory using the formalism of quasi-monomials and prove the usefulness of the 
method for the solutions of nontrivial difference equations. 
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CN ' In this note we consider the so called factorial polynomials [I], [2] by framing the relevant 

theory within the context of quasi-monomials [3]. The method we propose allows a significant 

simplification of the study of their properties and of their applications in the theory of multi-loop 

f"^ ■ Feynman integrals [4J. 

CN . The factorial polynomials (f.p.) 

C^ ■ ^"(") = r(x + i-n) (') 



< 



are quasi-monomials in the sense that they behave as ordinary monomials under the action of the 
operators 
C^ i M = xe-^" P = e^"-1, (2) 

nameljo 

M ipn (x) = ipn+1 (x) P ipn (x) = Tl (fn-lix) ■ (3) 

These operators are called multiplicative and derivative operators, respectively, and satisfy the 
commutation relation 



% 



1. (4) 
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Tlj- ■ The polynomials fn{x) satisfy the "differential" equation 

^: M P ipnix) = mpnix) , (5) 

that, according to eq. ([2]), is equivalent to the difference equation 

X [(pn{x) - ifnix - 1)] = rnp„{x) . (6) 

By defining the operator 

oo 

(t>{M) = Y,'^nM^, (7) 



n=0 

the identity 

M"¥Po(x) = M"l = 99„(x) (8) 

allows us to define a f.p.-based function as follows 

oo 

/(x) = </.(M)l = ^an<^„(x). (9) 

n=0 

For example, to the operator exp(AM) it is associated the f.p.-based function 

°° \n 
e{x, A) = e^*^ 1 = V — Mx) = (1 + A)^ (10) 

n=0 



Though not in the explicit realization of Eq. ((5)), the operators M and P have been introduced in ref. [3]- 
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that plays the same role of exponential function for ordinary monomials since it is an eigenfunction 
of the operator P with eigenvalue A 

Pe{x,X) = Xe{x,X)- (H) 

In terms of it, the f.p. -based cosine and sine functions can be defined as follows 

c{x) = - [e(x, z) + e(x, -i)] = J^ ^ M^'^ 1 = ^ i^ ^,,{x) 

s(x) = — [e{x, i) - e{x, -i)] = V , \~ ^ ,, M^^+^ 1 = V , \ ' ,, (/^2fe+i(x) . (12) 

fe=o ^ ' k=o ^ ' 



The f.p. -based function defined in eq. (jlOp is involved in the solution of the finite difference heat 
equation 

PrW{x,T) = dlW(x,T) W{x,^) = w{x) (13) 

where 

PrW{x,T) = W{x,T + \)-W(x,T). (14) 

The solution formally writeqj 

W{x,T) = &y.Y>{Mrdl)w[x)\r, (15) 

that, taking into account eq. ([TO]) , allows us to get the solution of eq. ([13]) in the form 

W{x,T) = e{T,dl)w(x). (16) 

The expansion in terms of factorial series is 

W{x^t) = Y^'^^w^^-\x), (17) 

n=0 

where the upper index (2n) denotes the order of the derivative and t^ = mTQ is the discrete time, 
corresponding to the definition given in eq. ()14p of the time derivative (tq = 1 because we have 
chosen a unity time interval). The sum run up to ?n, as a consequence of the fact that ^m+i {tti) = 0. 
The method could be usefully exploited to deal with discrete time Schrodinger equations [5j. 



In the same way, starting with the operators 

fcr("n-hA:)!2"+2fc 



(3n{M) = Y .. A. ^^^ M''^'^^, (18) 



k=0 



the f.p. -based Bessel functions can be defined as Bn{x) = /3(M) 1. In analogy with the ordinary 
cylindrical Bessel functions, the operators (fT8]l satisfy the recurrences 



2P/3„(M) = [/3„_i(M) -/3„+i(M) 

2n/3„(M) = M k_i(M) + /3„+i(M) . (19) 



that, in terms of finite difference, gives 



B^ix + 1) = B^x) + i [B.n-i{x) - B„+i(x)] 
2nBn{x) = x [Bn^iix - 1) + Bn+i{x - 1)] . (20) 



The sufHx r indicates that we consider only the space of this variable. 



By denoting with A^ a number operator such that N Bn{x) = nBn{x), the combination of eqs. I\W\} 
ahows the introduction of the following index shifting operators 

E+ = -P + M-^N, E^ = P + ]\/r'^N, (21) 

whose action on the operators in eq. (jlSp is given by 

i±/3„(M) = /3„±i(M) (22) 

from which it is easy to show that 

{MPf + (Af2 - n^)] /3n{M) = . (23) 

The corresponding difference equation writes 

2x{x-l)Bnix-2) -x{2x-l)Bn{x-l) + {x^ -n'^)Bn{x) = 0. (24) 



So far a differential equation has been translated into a corresponding difference equation using 
the realization ([2]) of the multiplicative and derivative operators. The correspondence between the 
solutions (i.e. their isospectrality) is ensured by the fact that the operators P, M, and 1 realize 
a Weyl algebra. Let us now consider the same problem from a reversed perspective, i.e., given a 
difference equation we obtain from it a corresponding differential equation in terms of derivative 
and multiplicative operators. 

Let us consider the following difference equation 

{ix + 2) F{x + 1) + 4x F{x - 1) - i8x-l)F{x) =0. (25) 

By performing the following substitutions 

e^^ ^ 1 + P , X ^ M (1 + P) , (26) 

we find 

^4Mp2 + 2 P + l) $(M) = (27) 



where <l>(Af)) 1 = F{x). The solution of the previous differential equation can be found by using 
the Frobenius method. Hence, we look for solutions of the type 

oo 

$(M) = ^afcM'=+^ (28) 

fe=0 

and the associated indicial equation is 

2aoc(2c-l) =0. (29) 

By using the solution c = of this equation, we get 

P(x) = c^(M) 1 = E fe(r ^' = E fe(r ^^(-) ■ (30) 

fc=0 ^ '' k=0 ^ '' 

The other solution, c = 1/2, implies the definition of the operator M^", whose meaning will be 
discussed in a forthcoming paper, where the problems associated with the fractional powers of 
operators of the type ^ will be considered. 

The method we have envisaged can also be extended to the non-homogeneous case. For example, 
the equation 

ay{x + l)+hy{x)= g{x) (31) 

can be cast in the form 

Q^{M) = g{M{l + P)) (32) 



where Q = aP + {a + b) and y{x) = ^(M) 1. The particular solution of this equation formally 
writes 

V'p(M) = Q-i<7(M(l + P)) (33) 

and the explicit form can be written as 

^p{M)= dse-'^g{M{l + P)) (34) 

Jo 

from which, taking into account that 

e"^g(M(l + P))=5((M + a)(l + P)), (35) 

we obtain 

poo 

jpp{M)= / ds e"" ("+'') 5((M- a s)(l + P)). (36) 

Jo 

As an example, in the case g{x) = x^, by using eqs. ([3]) and ([8]), one obtains: 

yp{x) = ^{{x-z)^-z{l-z)} (^z = -^y (37) 

This note has just been aimed at fixing the general rule of the use of factorial series in the theory 
of difference equations. More specific problems regarding applications like those treated in Ref. [4] 
will be discussed elsewhere. 
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